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Abstract. We review an inflationary scenario with the anisotropic expansion rate. 
An anisotropic inflationary universe can be realized by a vector field coupled with 
an inflaton, which can be regarded as a counter example to the cosmic no-hair 
conjecture. We show generality of anisotropic inflation and derive a universal property. 
We formulate cosmological perturbation theory in anisotropic inflation. Using the 
formalism, we show anisotropic inflation gives rise to the statistical anisotropy in 
O | primordial fluctuations. We also explain a method to test anisotropic inflation using 

the cosmic microwave background radiation (CMB). 
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1. Introduction 

It is well known that inflation elegantly solves the horizon and flatness problems. 
Moreover, inflation accounts for the origin of the large scale structure of the universe. 
The point in an inflationary scenario is that the exponential expansion of the universe 
erases any classical memory, which is often referred to as the cosmic no-hair conjecture. 
Because of this feature, quantum fluctuations are responsible for the origin of the large 
scale structure of the universe. Remarkably, the nature of the primordial fluctuations is 
understood by symmetry in inflation. 

• homogeneity 

First of all, we have to assume the initial homogeneity. Indeed, inflation does not 
commence with strong inhomogeneous initial conditions ( There is a possibility that 
this symmetry breaks down [U EJ [3].). 

• shift symmetry 

In order to have slow-roll inflation, we need a sufficiently flat potential. Hence, we 
have a shift symmetry in field space, <p(x) — > <p(x) + 0. Here, is a constant. 

• temporal de Sitter symmetry 

The metric for de Sitter spacetime reads 



where t is a constant. 

• spatial de Sitter symmetry 

Once inflation occurs, the cosmic no-hair conjecture tells us that the universe will 
be isotropized in a few Hubble expansion. 

Here, we assumed single field inflation with a standard kinetic term. 

The above symmetry determines the nature of primordial fluctuations. In general, 
we need n-point correlation functions to characterize the statistical nature of primordial 
fluctuations. However, the shift symmetry in field space implies suppression of 
nonlinearity and hence the Gaussian statistics of fluctuations. Thus, we need only 
2-point functions. As an example, we take curvature perturbations (. In a Fourier 
space, we have the power spectrum 



ds 2 = -alt 2 + e 2Ht (dx 2 + aly 2 + dz 2 ) , 
where the Hubble parameter H is constant. It is easy to find isometry 




t ->• t + t , x i ->■ e~ m x i , 



x l = {x,y,z} , 



<C(k 1 ),C(k 2 ) >=P(k 1; k 2 ), 



where k is the wavenumber vector. Moreover, the homogeneity constrains the power 
spectrum to be 

<C(ki),C(k 2 ) >=5(k 1 + k 2 )P(k 1 ). 
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Here, the delta function stems from the "momentum" conservation. The spatial de 
Sitter symmetry further constrains the power spectrum as 

< C(kO ,C(k 2 ) >= 5(k a + k 2 )P(h = Ikxl). 

Namely, the direction dependence is forbidden by the rotational symmetry. Finaly, the 
temporal de Sitter symmetry yields a scale invariant power spectrum 

P(k) = const. 

because of the scale invariance of spatial coordinates. These predictions are robust and 
universal in inflationary scenarios. In fact, the above predictions have been confirmed 
by CMB observations [4]. 

However, precision cosmology forces us to look at fine structure of the primordial 
fluctuations. In fact, since the universe is not exactly de Sitter, there exists violation of 
the temporal de Sitter symmetry, which leads to a slight tilt of the power spectrum. As 
the deviation from de Sitter expansion can be characterized by the slow roll parameter, 
the tilt should be of the order of the slow roll parameter. Similarly, we have violation 
of the shift symmetry, which leads to non-Gaussianity characterized by the slow roll 
parameter [5]. Along this line of thought, it is natural to expect violation of the spatial 
de Sitter symmetry, which would lead to the statistical anisotropy. 

From the observational point of view, a lot of anomalies indicating the statistical 
anisotropy are reported although its statistical significance is uncertain ( see [6] and 
references therein). Motivated by those observations, there are many theoretical 
proposals to realize the statistical anisotropy [7J EJ |9l [lOj [HJ [12l ]. In these works, 
however, a consistent theoretical framework including the backreaction of vector fields 
seems to be obscure. As another line of research, there exist challenges to the cosmic 
no-hair conjecture |23j EH ESI EHJ E3 ]• If we can evade the cosmic no-hair conjecture, 
we would have the statistical anisotropy. Unfortunately, it turns out that these models 
suffer from either the instability, or a fine tuning problem, or a naturalness problem [32J . 
Recently, however, stable anisotropic inflation models are found in the context of 
supergravity, which gives rise to a counter example to the cosmic no-hair conjecture [33J. 

In this review, we explain how anisotropic inflation can be realized in supergravity 
and show it implies violation of spatial de Sitter symmetry and hence leads to the 
statistical anisotropy in primordial fluctuations. Recall that the bosonic sector of the 
supergravity action is given by 

S = J d 4 i^ - GijdPftdurft - e K2R (c^DiWDjW - 3k 2 WW) 



where Gjj = dK/dftdft, DiW = dW/dft + tf{dK/d<jf)W , K(<p,(f>) and W{<j>) are the 
Kaler potential and the super potential, respectively. There is also a kinetic term for 
gauge fields with gauge kinetic functions f a b- The cosmological role of K and W are well 
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discussed so far. However, the gauge kinetic function has been overlooked in cosmology, 
at least in the context of inflationary scenarios. In this review, we clarify the role of the 
gauge kinetic functions in inflation. 

Let us summarize main results here. First of all, we find anisotropic inflation is an 
attractor in supergravity with a wide rage of gauge kinetic functions. The metric during 
inflation approximately reads 

ds 2 = -alt 2 + e 2Ht [e-^alx 2 + e m (aly 2 + alz 2 )] , (1.3) 

where H and E describe the average expansion rate and the anisotropic expansion rate, 
respectively. The degree of the anisotropy H/H is universally given by the following 
formula 

S 1 H 

H = 3 l€H > eH = -JP> (L4) 
where I is the model parameter taking values < I < 1. The point is that the degree 
of anisotropy at most of the order of the slow-roll parameter e#. In this scenario, we 
have the statistical anisotropy of the form [31] 

P(k) = P(k) [l + ^sin 2 fl] , (1.5) 

where P(k) is the isotropic part of the power spectrum -P(k) and 9 is the angle between 
the preferred direction and the wavenumber vector of fluctuations. The amplitude of 
anisotropy g* can be calculated using the standard perturbation theory. The anisotropy 
in curvature perturbations is given by 

g s = 2AIN 2 (k) (1.6) 

and that in gravitational waves reads 

g t = QIe H N 2 (k) . (1.7) 

Here, / is a model parameter appeared in the anisotropy formula fll.4jl and N(k) is 
the e-folding number from the horizon exit of fluctuations to the end of inflation. 
Remarkably, there exists a difference between the anisotropy in curvature and tensor 
perturbations. There is also the cross-correlation between the curvature perturbations 
C and gravitational waves h given by 

r c = = -2AV2Ie H N 2 (k). (1.8) 

We find the consistency relations between these observables 

% = t H g s , r c = -V2e H g s . (1.9) 

This allows us to give a model independent test of anisotropic inflation. Of course, we 
can use each observables to constrain gauge kinetic functions. Indeed, we give the first 
cosmological constraint on the gauge kinetic function: 

where we used the result in [35]. Note that / is determined once the gauge kinetic 
functions are given. 
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The organization of the paper is as follows. In section II, we show power-law 
inflation is not necessarily an attractor in the presence of gauge kinetic function. Instead, 
we show anisotropic inflation could become an attractor for a wide range of parameters 
in models. In section III, we argue the generality of the model and derive a universal 
relation. In section IV, we develop cosmological perturbation theory in anisotropic 
inflation and calculate the statistical anisotropy in primordial fluctuations. In section 
V, we discuss observational test of anisotropic inflation using the CMB. The final section 
is devoted to summary and future prospects. 

2. Primordial Magnetic Fields, Backreaction, and Anisotropic Inflation 

In this section, we first recall a standard mechanism for generating primordial magnetic 
fields during inflation. There, a vector field coupled to an inflaton is introduced. Here, 
we consider a simple model with exponential potential and gauge kinetic functions. Then 
we obtain exact power-law inflation. It is possible to generate primordial magnetic fields 
in this set up. For many cases, the backreaction of the vector field cannnot be negligible. 
Indeed, it turns out that the backreaction yields a new type of cosmological solutions, 
namely, exact anisotropic power-law inflation. 

2.1. Power-law Inflation and Primordial Magnetic Fields 
Let us start with the action 



where k 2 is the reduced gravitational constant, g is the determinant of the metric, R is 
the Ricci scalar, V(<f>) is a potential for an inflaton <fi. In order to find exact solutions, 
we take the exponential potential 



where e a is the scale factor. Let us seek isotropic power-law solutions by putting the 
ansatz 





V{4>) = V e x ^ , 

where V and A are parameters. It is natural to take the isotropic metric 



(2.2) 



fa* = _ df 2 + e 2a(t) [ dx 2 + dy 2 + ^2 



(2.3) 



« = Clogt , 



= f log t + <f) • 



(2.4) 



Then, we obtain the solutions 

A 2 ' ? A 
In this case, we have the spacetime 




(2.5) 




(2.6) 



Thus, for A <C 1, we have power-law inflation. 
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Now, we consider primordial magnetic fields in this inflationary background. We 
introduce a vector field A^ whose kinetic term is coupled to the inflaton field <fi 



S 



d 4 x\/—g 



where F^ v is the field strength of the vector field defined by F, 



d„A„ 



(2.7) 
dyA^, and 

f(4>) is a coupling function of the vector field. We emphasize that this kind of model is 
quite natural in the context of the supergravity [36]. Many years ago, Ratra considered 
the exponential gauge kinetic function [37] 

70) = (2.8) 

and concluded that the primordial magnetic fields can be generated due to the gauge 
kinetic function. The result implies there exists a contribution of the vector field 
to the energy density in the universe. In the subsequent two subsections, we would 
like to show the backreaction of the vector field leads to the anisotropic inflationary 
power-law solutions and prove that anisotropic inflation is actually an attractor [38J. 
Thus, the models producing primordial magnetic fields naturally lead to anisotropic 
inflation [351 |4"0"] . Here, we should emphasize that the existence of anisotropic inflation 
attractor shows that the backreaction from the vector field does not necessarily destroy 
inflation |39|, as is often assumed in literature [41]. 



2.2. Backreaction and Anisotropic Power-law Inflation 

Now, we take into account the backreaction of the vector field. We can expect interesting 
effects due to the backreaction because the coupling of the vector field to the inflaton 
produces an effective potential for the inflaton [121 S3]- Indeed, we will see that there 
exist exact anisotropic solutions [38] on the contrary to the expectation from the cosmic 
no-hair theorem [M 1431 I4U1 147] . 

The action we should consider is given by 



S 



d 4 Xy/—g 



1 



1 



V(<P) - -f^fF^ 



,(2.9) 



where V((j)) and f((f>) are given by (12. 2p and (12. 8p . respectively. Without loosing the 
generality, one can take x-axis in the direction of the vector field. Using the gauge 
invar iance, we can express the vector field as 

= v(t)dx . (2.10) 

Thus, there exists the rotational symmetry in the y-z plane. Given this configuration, 
it is convenient to parameterize the metric as follows: 



ds 2 



M{t) 2 dt 2 + e 2a ^ \e-^dx 2 + e 2 ^\e 2V ~ 3 ^dy 2 + e - 2V ~ 3 ^dz 2 )} , (2.11 



where e a , a and cr_ are an isotropic scale factor and spatial shears, respectively. Here, 
the lapse function M is introduced to obtain the Hamiltonian constraint. With the 
above ansatz, the action becomes 



S 



d 4 xj^e 3a 



1 



— {-a 2 + a 2 + &l) + - 



> 2 - U 2 V{<t>) + i/(0)V e - 2a W +4CT « 



(2.12) 
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where an overdot denotes a derivative with respect to the physical time t. First, its 
variation with respect to o~_ yields 

o_ = -3a&- . (2.13) 

This gives o\_ oc e _3a , hence, the anisotropy in the y-z plane rapidly decays as the 
universe expands. Hereafter, for simplicity, we assume cr_ = and set the metric to be 



als' 



-alt 2 + e 2a ^ 



-4(7 



x + e 



2a(t) 



(aly 2 + alz 2 ) 



(2.14) 



Next, the equation of motion for v is easily solved as 

v = f(cf>)- 2 e- a -^ PA , (2.15) 

where pa is a constant of integration. Taking the variation of the action with respect to 
A/", a, cr and <p an d substituting the solution fl 2 . 1 5 [) into them, we obtain the following 
basic equations: 



a 



a 



a 



<t z + 



l^ 2 + vw + ^n<pr 2 e-^ 



- 3d 2 + K 2 V(<f>) + ^i/(0)- 2 e - 4Q - 4CT 



6 



3a& 



K 2 p 2 A 



-2 e -4o-4<7 



(2.16) 

(2.17) 

(2.18) 
(2.19) 



= -3d0-^+pi/(0)-%e- 

where the subscript in denotes a derivative with respect to 0. Let us check whether 
inflation occurs in this model. Using Eqs. (I2.16P and (I2.17p . the equation for acceleration 
of the cosmic expansion is given by 



a + cY 



Z<7 — 



K 



2 

+ 0L y _PA f -2 -4a-4a ^ ^ 

e* 3 T 3 [ 2 J K J 

We see that the potential energy of the inflaton needs to be dominant and the energy 
density of the vector field 
1 



Pv = -^P A f{4>)' ~< 



(2.21) 



and the shear £ = a should be subdominant for inflation to occur. 
In order to find exact solutions, we take the power-law ansatz 

a = (logt, a = r]\ogt, K(f) = flog* + 0o . (2.22) 

From the hamiltonian constraint equation (12.161) . we get two relations 

Af = -2 , pi + 2C + 2 V = 1 (2.23) 

to have the same time dependence for each term. The latter relation is necessary only 
in the non-trivial vector case, p A ^ 0. Then, for the amplitudes to be balanced, we need 

, (2.24) 

o u 

where we have defined variables 



C 2 + r? + -£ 2 + -u + -w 
S ' 6^ 3 6 



u 



K 2 V Q e> 



w 



lie- 



-2p<p 



(2.25) 
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The equation for the scale factor (12.171) under Eq. ( I2.23P yields 

-( + 3( 2 -u--w = . (2.26) 
6 

Similarly, the equation for the anisotropy ( I2.18P gives 

-V + Kv-^w = . (2.27) 
Finally, from the equation for the scalar (" 12 . 19j) . we obtain 

- f + 3C£ + \u - pw = . (2.28) 
Using Eqs. ( I2.23p . (I2.26P and (I2.27p . we can solve u and w as 

^ 2 -i<-^+; + ^ <"■> 

Substituting these results into Eq. (I2.28p . we obtain 

(3C - 1) [6A (A + 2p) C - (A 2 + 8pA + 12p 2 + 8)] = . (2.31) 

In the case of ( = 1/3, we have u = w = 0. Hence, it is not our desired solution. Thus, 
we have to choose 

_ A 2 + 8pA + 12p 2 + 8 
C " 6A(A + 2p) • [Z - 6Z) 

Substituting this result into Eq. ( I2.26p . we obtain 

» + 2pX-* (2 . 33) 

' 3A(A + 2p) v ' 

This clearly shows the existence of the anisotropy in the expansion. From Eq. ( I2.23p . 
we have 

e=~ ■ (2.34) 

Finally, Eqs. (12T29|) and (12~30|) reduce to 

(pA + 2p 2 + 2)(-A 2 + 4pA + 12p 2 + 8) 

2A 2 (A + 2p) 2 1 ' ' 

and 

2A 2 (A + 2p) 2 v ; 

Note that Eq. (12. 24ft is automatically satisfied. Thus, we have obtained anisotropic 
power-law solutions. 

Recalling the definition (I2.22p . we see C, ^> 1 is necessary for inflation. From the 
solution (I2.32p . it turns out that this requirement can be achieved by assuming A p. 
For these cases, u is always positive. Since w should be also positive, we have the 
condition 

A 2 + 2pA > 4 . (2.37) 
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Hence, p must be much larger than one. Now, the spacetime reads 

ds i = - d t 2 + t 2 ^dx 2 + t 2 ^ (dy 2 + alz 2 ) . (2.38) 

The average expansion rate is determined by ( and the average slow roll parameter is 
given by 

tH ~ H 2 ~ X 2 + 8p\ + 12p 2 + 8 ' { ' 

where we have defined H = a. In the limit A 1 and /)> 1, this reduces to en = A/p. 
Now, the anisotropy is characterized by 

£ _ a 2(A 2 + 2pA - 4) 



H a A 2 + 8pA + 12p 2 + 8 ' 
From Eq. f!2.39j) and (12.401) . we obtain a relation 



(2.40) 



El T A 2 + 2pA-4 . A . 



It is possible to write / as 

_ c - 1 A 2 + 2pA 



(2.42) 



c ' 4 

Then, it is apparent / takes a value in the range < / < 1. We see the anisotropy is 
positive and proportional to the slow roll parameter en- 

Although the anisotropy is always small, it persists during inflation. Clearly these 
exact solutions give rise to counter examples to the cosmic no-hair conjecture. We 
should note that the cosmological constant is assumed in the cosmic no-hair theorem 
presented by Wald [44J . In the case of isolated vacuum energy, the inflaton can mimic the 
cosmological constant. However, in the presence of the non-trivial coupling between the 
inflaton and the vector field, the cosmic no-hair theorem cannot be applicable anymore. 

2.3. Anisotropic Inflation as an Attractor 

In the previous subsections, we found both isotropic and anisotropic power-law solutions 
exist. In this subsection, we will investigate the phase space structure. Then, we will 
see which one is dynamically selected. 

Let us use e-folding number as a time coordinate da = adt. It is convenient to 
define dimensionless variables 

X = - , Y = , Z = K fU)e- a+2a - . (2.43) 
a a a 

With these definitions, we can write the hamiltonian constraint equation as 

- K 2 Y- = 3(X 2 - 1) + -Y 2 + -Z 2 . (2.44) 
a 2 v ; 2 2 v ; 

Since we are considering a positive potential, we have the inequality 

3(X 2 - 1) + W 2 + ^Z 2 < . (2.45) 
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Using the hamiltonian constraint (12.44j) . we can eliminate from the equations of 
motion. Thus, the equations of motion can be reduced to the autonomous form: 



^ = iz 2 (X + l)+x{3(X 2 -l) + ir 2 } 
(F + A)(3(X 2 -1) 



dY 
da 



-Y Z 



P + 



— = Z 

da 



3 (X s 



-Y — pY + 1 — 2X + -Z z 



(2.46) 
(2.47) 
(2.48) 



Therefore, we have a 3-dimensional space with a constraint (I2.45p . A fixed point in this 
phase space is defined by dX/da = dY/da = dZ/da = 0. 

First, we seek the isotropic fixed point X = 0. From Eq. (12.461) . we see Z = 0. 
The remaining equation (12.471) yields Y = —A or Y 2 = 6. The latter solution does not 
satisfy the constraint (12.451) . Thus, the isotropic fixed point becomes 

(X,F,Z) = (0,-A,0). (2.49) 

This fixed point corresponds to the isotropic power-law solution (12. 61) . Indeed, one can 
check that the solution (12. 5p leads to the above fixed point. 

Apparently, Z = and 6X 2 + Y 2 = 6 give a fixed curve. However, this contradicts 
the constraint (I2.45j) . 

Now, let us find an anisotropic fixed point. From Eqs. f!2.46j) and (I2.47p . we have 

A^ 



Y 



Eq. (PX461) gives 



3p4 



3X 
X + l 



X- A . 



3(X 2 - 1) + -Y 2 



Using the above results in Eq. (I2.48[) . we have 



X-2 A 2 + 2pA-4 



. 



(2.50) 

(2.51) 
(2.52) 



(X - 2) [(A 2 + 8pA + 12p 2 4 

The solution X = 2 does not make sense because it implies Z 2 = —18 — 36p 2 < by 
Eqs. (I2.50p and (I2.5ip . Thus, an anisotropic fixed point is expressed by 

2 (A 2 4 2pA - 4) 



X — 

A 2 4 8pA 4 12p 2 4 8 ' 
Substituting this result into Eq. (I2.50p . we obtain 

12 (A 4 2p) 



Y 

Eq. fl2T5T]) yields 
Z 2 



A 2 4 8pA 4 12p 2 4 8 ' 
18 (A 2 4 2pA - 4) (-A 2 4 4pA 4 12p 2 4 



(2.53) 



(2.54) 



(2.55) 



(A 2 4 8pA 4 12p 2 4 8) 2 

Note that from the last equation, we find that A 2 4 2pA > 4 is required for this fixed 
point to exist under inflation A p. It is not difficult to confirm that this fixed point 
corresponds to the anisotropic power-law solution (12.381) . 
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Next, we examine the linear stability of the fixed points. The linearized equations 
for Eqs. fl2T46|) . ( HAT} . f|2T48|) are given by 

= (\z 2 + 9X 2 + \y 2 -Z)8X + XY5Y + \ {X + 1) Z5Z , (2.56) 
da \6 I J 3 

dSY „ f ,„„ \ 1 „ 1 



da 



6X (y + A) <5X + 1 3 (X 2 - l) + -Y 2 + y (Y + A) + -Z 2 } <5F 
+ ( \y + 2p + \) Z5Z , (2.57) 



dSZ 
da 



.3 

2(3X - l)Z8X + (Y-p) Z5Y 

+ (^3X 2 + iy 2 + Z 2 - 2X - pY - 2^ SZ . (2.58) 
In the case of the isotropic fixed point Eq. (12.491) , these equations reduce to 

llSX 'h*-3)6X, (2.59) 



da 

£-(5*-»)". <-) 

d5Z 



1a 2 -2 + p A 



.. _ , . SZ. (2.61) 

da 12 

We see that the coefficient in the right hand side of above equations becomes negative 
when A 2 + 2pA < 4 during inflation A < 1, which means the isotropic fixed point is 
an attractor under these conditions and the isotropic fixed point becomes stable in this 
parameter region. In the opposite case, A 2 + 2pA > 4, the fixed point becomes a saddle 
point and unstable. In the latter cases, if there exist the vector field whatever small it 
is, the vector field destabilize isotropic inflation. 

Now we are interested in the fate of trajectories around the unstable isotropic fixed 
point. We will see that those trajectories converge to an anisotropic fixed point. Since 
we are considering the inflationary universe A<1, the condition A 2 + 2pA > 4 implies 
p ^> 1. Under these conditions, we can approximately write down the linear equations 
as 

dSX 



da 
d5Y 



35X, (2.62) 



da 

d5Z I 



3<5Y + a/6(A 2 + 2pA - A)5Z , (2.63) 



-j6(X 2 + 2p\-A)5Y . (2.64) 



da 2 
The stability can be analyzed by setting 

SX = e wa 5X , 5Y = e ua 5Y , 5Z = e wa 5Z . (2.65) 

Then we find the eigenvalues u are given by 



w = -3 , -| ± ^3(A2 + 2pA - 4) - H . (2.66) 

As the eigenvalues have negative real part, the anisotropic fixed point is stable. Thus, 
the end point of trajectories around the unstable isotropic power-law inflation must be 
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Z 




isotropic anisotropy 
inflation 

Figure 1. The phase flow in X-Y-Z space is shown for A = 0.1, p = 50. The 
trajectories converge to the anisotropic fixed point. 



anisotropic power-law inflation. In FiglTJ we depicted the phase flow in X-Y-Z space 
for A = 0.1, p = 50. We see that the trajectories converge to the anisotropic fixed 
point indicated by yellow circle. The isotropic fixed point indicated by orange circle is 
a saddle point which is an attractor only on Z = plane. Thus, anisotropic power-law 
inflation is an attractor solution for parameters satisfying A 2 + 2pA > 4 [38J. 



3. Anisotropic Inflation: Generality and Universality 

3.1. Generality 

Next we want to clarify the generality of anisotropic inflation. First, we need to look at 
the ratio of the shear to the expansion rate S/i? to characterize the anisotropy of the 
inflationary universe. Notice that Eq.f l2.18]) reads 

2k 2 

E = -3#E + — p v . (3.1) 
o 

If the anisotropy converges to a value, i.e. £ becomes negligible, the terminal value 
should be given by 

§ = i*- K = mr (3 - 2) 

where we used the slow roll equation 

H 2 = |W) , (3-3) 
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which is derived from Eg. (12.16j) . 

In order to realize the above situation, p v must be almost constant. Assuming that 
the vector field is subdominant in the evolution equation of the inflaton field Eq.f l2.19l) 
and conventional single field slow-roll inflation is realized, one can show the coupling 
function /(</>) should be proportional to e~ 2a to keep p v almost constant. In the slow 
roll phase, e-folding number a is related to the inflaton field (ft as da = —^V^dc^/V^ 
as usual. Then, the functional form of /(</>) is determined as 

/ (0)= e - 2a = e 2K2/ ^. (3.4) 

For the polynomial potential V oc n , for example, we have / = e~^~ . 

The above case is, in a sense, a critical one. What we want to consider is super- 
critical cases. For simplicity, we parameterize /(</>) by 

M=e*'*J%» (3.5) 

where c is a constant parameter. Now, we look at what happens when c > 1. Note that 
Eq. (13.51) can be written as 

4 = ■ (3.6) 

Then, the condition c > 1 can be promoted to the condition 

> 1 ' (a7) 

Thus, any functional pairs / and V which satisfies (13.71) in some range could produce 
the vector-hair during inflation. The equation for the inflaton becomes 

2c 



-K 



(3.8) 



where we have defined the slow-roll parameter 

In this case, if the vector field is initially small 1Z <C ey/2c, then the conventional single 
field slow-roll inflation is realized. During this stage / oc e~ 2ca and the vector field 
grows as p v oc e 4 ( c-1 ) a . Therefore, the vector field eventually becomes relevant to the 
inflaton dynamics Eg. (13 . 8 p . Nevertheless, the accelerating expansion of the universe 
will continue. The point is that TZ cannot exceed ey/2c. In fact, if TZ exceeds ey/2c, 
the inflaton field <fi does not roll down, which makes p v = p^/(0)~ 2 e _4a_4<7 /2 decrease. 
Hence, p v V(<f>) always holds. In this way, there appears an attractor where the 
inflation continues even when the vector field affects the inflaton dynamics. 
The inflaton dynamics is determined by solving the slow-roll equation: 

- 3d0 -Vf + p 2 A r 3 Ue- 4a - Aa = . (3.10) 

Using the slow- roll equation (13. 3ft . this yields 

# _ _ ^V^_ + 2c P 2 ie -4a-4<x-W/^# _ (3 n) 



da a k 2 V V« 
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This can be integrated by neglecting the evolutions of V, V^, a as 



Aa+4a+Ac K 2 / _ 2c 2 p 2 A K 2 V 



1 + fie- 4(c - 1)Q+4<J l , (3.12) 



'c-lV 2 

where Q is a constant of integration. Substituting this into the slow-roll equation 
Eq. (13. Ill) . we obtain 



V., c- I V„ 



1 + fie- 4(c - 1)Q+4CT l 1 . (3.13) 



da n 2 V c k 2 V 
Initially a — > — oo, the second term can be neglected. While, in the future a — > oo, the 
term containing Q disappears. This clearly shows a transition from the conventional 
single field slow-roll inflationary phase, where 

holds, to what we refer to as the second inflationary phase, where the vector field is 
relevant to the inflaton dynamics and the inflaton gets 1/c times slower as 

da c k 1 V 

In the second inflationary phase, we can use the formula (13.121) dicarding Q term and 
rewrite the energy density of the vector field as 

p\ —Aa—Au—AcK 2, [ T7-d(f> lc— 1 , N , , „ N 

Pv = Y e * = 2^ eyV ^ ' (3 - 16) 

which yields the anisotropy 

£ _ 2 lc- 1 
H~3 3~ 2 

Moreover, from Eqs. (l2.16l) and (12.171) . the slow- roll parameter defined in terms of the 
scale factor becomes 

a 1/1 2; , 2 2 \1 
*h = -# = -~ 2 {- r <P ~ r P«) = ~ey , (3.18) 

where we neglected the anisotropy and used relations ( 13.31) and (13.151) . Thus we have a 
remarkable result [33] 

E lc- 1 , 

77 = 3—— (3 ' 19) 

Therefore, for a broad class of potential and gauge kinetic functions, there exist 
anisotropic inflationary solutions. 



17 = o n = o— —ev ■ (3.17) 



3.2. Example: Chaotic Inflation 

In order to make the statement concrete, we consider chaotic inflation with the potential 

V&) = \m 2 <i> 2 , (3.20) 

where m is the mass of the inflaton. For this potential, the coupling function becomes 
/(0) = e^ 2 * 2 / 2 . It is instructive to see what happens by solving Eqs. (12.16j) - (12. 19j) 
numerically [33J. 



Statistical Anisotropy from Anisotropic Inflation 



15 




c=z 



anisotropic inflation 



isotropic 
inflation 



-19-005 ' 1 1 1 ' 1 1 1 

-2 2 4 S 8 10 12 

4> 



Figure 2. Phase flow for <j> is shown. Here, we took the parameters c = 2 and 
urn = 10~ 5 . We also took initial conditions fa — 12 and ^ = 0. There are two 
different slow-roll phases. The transition occurs around K<p — 9. 




Figure 3. Evolutions of the anisotropy Yi/H for various c with respect to the e- folding 
number are shown. One can see the attractor behavior of the anisotropy. 



In Fig. [21 we have shown the phase flow in — space where we can see two 
slow-roll phases. The first one is the conventional inflationary phase and the second one 
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is the anisotropic inflationary phase. As usual, inflation ends with oscillation around the 
bottom of the potential. This tells us that isotropic inflation corresponds to a saddle 
point and ansiotropic inflation is an attractor in the slow roll phase. In contrast to 
anisotropic power-law inflation, the above chaotic anisotropic inflation is a transient 
phase in the whole phase space. 

In Figj3l we have calculated the evolution of the anisotropy H/H = a /a for various 
parameters c under the initial conditions y/cK<pi = 17. As expected, all of solutions 
show a rapid growth of anisotropy in the first slow-roll phase which corresponds to 
the conventional inflation. However, the growth of the anisotropy eventually stops at 
the order of the slow roll parameter. Notice that this attractor like behavior is not so 
sensitive to the parameter c. As one can see there is a sufficient amount of e-folding 
number during an anisotropic inflation. 

3.3. Universality 

As we have seen, the anisotropy satisfies the inequality 



This inequality holds universally for any potential functions. This result is reasonable 
because the cosmic no-hair theorem holds for a gravity system with a positive 
cosmological constant dominating the universe. Since the deviation from the exact 
de Sitter is characterized by the slow roll parameter, the deviation from the isotropic 
expansion must be proportional to the slow roll parameter. 

3.4- Variety of Models 

There are many models which realize anisotropic inflation. We can generalize single field 
inflation to multi- field inflation models [HI H9]. Indeed, in almost all kind of models 
including small field and hybrid inflation, there exists anisotropic inflation. Actually, a 
more wide range of anisotropic inflationary models are discussed in [50]. We can extend 
the standard kinetic term to the Born-Infeld type j5TJ [52], [53]. In this direction, we may 
find a stringy realization of anisotropic inflation. We can introduce a mass term to the 
vector field, that is, a vector curvaton [511 [21] . It is interesting to study cosmological 
consequences of the vector curvaton scenario in detail. It is also possible to extend 
the model to non-abelian gauge fields [55]. In this case, we have more complicated 
dynamics which would lead to interesting phenomenology. Interestingly, there are other 
non-abelian gauge field models [561 [571 EE] ■ It is interesting to extend the analysis in this 
paper to other Bianchi type models [59]. We can consider inflation with multi- vector 
fields [60] and other tensor fields to realize anisotropic inflation. 





Statistical Anisotropy from Anisotropic Inflation 



17 



4. Statistical Anisotropy in Primordial Fluctuations 

So far, we have shown the existence of the anisotropic hair in a variety of inflationary 
models. It is interesting to see how to test anisotropic inflation using observations. 
To this end, we need the information of fluctuations in anisotropic inflation. Since 
the background is anisotropically expanding, we cannot use the standard cosmological 
perturbation theory [6TJ [62j [63j [MJ [65j EE E7J [68]. In this section, we classify 
perturbations under the 2-dimensional rotational symmetry and obtain the quadratic 
actions for 2-dimensional scalar and vector sectors. In order to grasp the meaning of 
variables, we start with the isotropic case and make a gauge transformation from the 
flat slicing gauge to the appropriate gauge for 2-dimensional classification. Then, the 
resultant gauge can be promoted to the anisotropic spacetime. The gauge we have 
chosen makes the analysis and the interpretation easier. Once the gauge is fixed, it is 
straightforward to calculate the quadratic action. 



4-1. Gauge Fixing and Classification of perturbations 



First, we start with the spatially homogeneous and isotropic universe. For simplicity, 
we consider flat space. 



ds 2 = a 2 (rj) —drf + 8ijdx l dx^ 



(4.1) 



where we took a conformal time rj. In that case, we can use 3-dimensional rotational 
symmetry to classify the perturbed metric. When we want to have the diagonal 
quadratic action, we take the following gauge 



ds' 



(4.2) 



-(1 + 2A)drj 2 + 2(£; + VAdrjdx* + (<% + h^d^dx^ , 

where we imposed V% } % = and hijj = ha = 0. If we ignore vector and tensor 
perturbations V i: hij, the above gauge is called the flat slicing gauge. Now, let us move 
on to the Fourier space. Since there exists 3-dimensional rotation symmetry, we can take 
a wavenumber vector to be k = (k, 0, 0). Then, the perturbed metric has the following 
components: 



8g. 



in/ 



( -2a 2 A a 2 B r 
* 



V 



a 2 V 2 



a 2 h. 



a 2 V 3 


—a 2 h > 
-a 2 h 



\ 



* is symmetric part. (4.3) 



+ / 



Here, we utilized the special choice k = (k, 0, 0) to solve the constraints = and 



j = ha = 0. With the same reason, only B x remains. We defined h yz 



-h h 



—h zz = h + . Now, we will pretend that we have only 2-dimensional rotation symmetry in 
y-z plane. In that case, at best, we can take k = (k x , k y , 0). Hence, we make a rotation 
in the x — y-plane so that the wavenumber vector becomes k = (k x , k y , 0). 



/ k x \ 



k i, 



( k x 


- ky 


\ 


( 


k \ 


ky 


k x 









V o 





) 


V 





(4.4) 
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where we have a relation k 2 
becomes 



k x + k 2 . Under this rotation, the perturbed metric 



/ -2a 2 A 



Sg 



/IV 



!fa 2 B }X - k fa 2 V 2 !fa 2 B, x + >fa 2 V 2 



a 2 V 3 



n 2 ^h 

a p 1 1- 



\ 



-a 2 ^h, 
a 2 §h+ 
* 



k 

*2 k x 



a k h x 



—a 2 h 



(4.5) 



+ / 



To simplify the perturbations, we can make use of gauge transformation 

Sg^, -> Sg^ + £ w „ + ^ , (4.6) 

where the semicolon denotes the covariant derivative with respect to the background 
metric. Taking the parameter 

A',: 







2ik 2 



h. 



2ik 2 



-——h 



we obtain 



/ -2a 2 A 



!fa 2 B tX + 
a 2 h. 



V 



!fa 2 B, x + 



a 2 h A 



a 2 V 3 + ■ 

Ky 



—a 2 h~ 



(4.7) 



where we have omitted some unimportant parts. It should be noted that we did not 
change slicing but performed only the spatial coordinate transformation. Therefore, we 
are still working in the flat slicing where the 3-dimensional scalar curvature vanishes. 

In our anisotropic inflation models, the available symmetry is actually small. The 
background metric is given by 

a(7]) 2 (-dr] 2 + alx 2 ) + W + dz 2 ) } (4.8) 



olsl 



that is, a = e a ~ 2cr ,b = e a+r7 ,dr] = alt/ a. Notice that the conformal time in anisotropic 
inflation is the conformal time in 2-dimensional part (r],x). Even in this anisotropic 
spacetime, as we have done in f)4.7p . one can take the following gauge 



\ 



8 goo 


Sgox 


tigoy 


Sgoz 


* 


Sg xx 





Sg xz 


* 


* 


5g yy 





* 


* 


* 


Sg zz 



(4.9) 



where we can impose further conditions so that the perturbed metric goes back to (I4.7I) 
in the isotropic limit. 

One can classify the perturbed metric using the rotational symmetry in y — z-plane. 
In 2-dimensional flat space, an arbitrary vector m a where a = y, z can be decomposed 
into the scalar part m a a ^ and the vector part m° = 0. Since there exists no tensor 
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part in 2-dimensions, 2- dimensional tensor can be constructed from the 2-dimensional 
vector. Because of the symmetry, the scalar and vector parts are not mixed in the 
equations. Thus, the metric perturbations can be classified into the scalar sector and 
the vector sector. Thanks to the symmetry in the y — z plane, without loss of generality, 
we can take the wavenumber vector to be k = (k x ,k y ,0). Hence, the vector sector in 
2-dimensional classification can be identified as 5go z , 5g xz in the above perturbed metric. 
The remaining components 5g 00 , 5g 0x , 5g 0y , 5g xx , 5g yy , 5g zz belong to the scalar sector. 



4-1.1. 2d vector sector Thus, the perturbations that belong to 2d vector perturbations, 
can be written as 

/ n n n b 2 /3 3 \ 

b 2 T 





\ 





* 
* 
* 






* 
* 







* 



(4.10) 



where we have incorporated the anisotropy while keeping the spatial scalar curvature 
to be zero. As to the vector field, we can take 



^vector 



(0,0,0,1?) 



(4.11) 



Note that we have no residual gauge transformation and, in particular, I? is a gauge 
invariant under abelian gauge transformations. And, as we have seen in (14. 7p . T 
corresponds to the cross-mode polarization of gravitational waves in the isotropic limit 
a = b. 

Using this gauge, we can calculate the quadratic action as 



S 



vector 



drjafx 



o 4 b 2 



o 4 f 2 v'b 2 



2a 



^2 

— r 2 

4 * 



b A f 2 a 2 
4a 2 26 2 ' y 



1 



f 2 f 2 v'b 2 



(4.12) 



2" > x 2 a' 

Since the perturbed shift function /3 3 does not have a time derivative, it is not dynamical. 
There are two physical degrees of freedom V and D in this 2-dimensional vector sector. 



4-1-2. 2d scalar sector For the 2-dimensional scalar sector, we define the metric 
perturbations 



scalar 



-2a 2 $ 



V 



* 
* 



a(3i 
2a 2 G 



* 
* 



a(3 2 


2b 2 G 



\ 



(4.13) 








-26 2 G J 

where we have kept the spatial scalar curvature vanishing. The scalar perturbation 
will be represented by 5(f). The variable G and 5(f) are the gauge invariant variables 



Statistical Anisotropy from Anisotropic Inflation 



20 



that correspond to the plus mode of gravitational waves and the scalar perturbations, 
respectively, in the isotropic limit a = b. And, we set the perturbed vector to be 

M scalar = (Mo)0)J)0) ? (414) 

where we have fixed the abelian gauge by putting the longitudinal component to be 
zero. From these ansatz, we can calculate the quadratic action as 



S 



scalar 



al 3 xdr] 



i2 r2 l2 

7^fSAl x + J -8Al y + -f 2 v' (G + $) 5A, X - f 2 J'5A 0! y 



b 2 11 bb' b 2 1 

2-ff (t> v'5cj ) 5A , x + - A (3l y - -fa^y + 2— $ ^ - - + 



+ a 



+ 



If 



V 



o I a 2 



~ \fJ% + b 2 G' 2 - a 2 G 2 y - b 2 G 2 x + h?^' 2 - pfi - p<j>% - \a 2 b 2 V H ^ 2 
+ f~ 2 (/J + ffn) v'V - a 2 b 2 V& + f- 2 f 2 v' 2 G 2 - 2a% 2 V$G - 2bb'$'G 



2a 2 ' 

(G + $) + 6 W (G - $) 



(4.15) 



Here, S' scalar consists of /3i, /3 2 ,G,5A ,5(J) and J. Among them, <£>,/3i,/3 2 and 5A are 
non-dynamical and can be eliminated. 

In order to calculate the statistical properties of primordial fluctuations from 
anisotropic inflation [69] [70] [71] [72] , we need to reduce the action to the one for physical 
variables. Then, we can quantize the system and specify the vacuum state. We analyze 
the vector sector and the scalar sector, separately. 



4-2. Action in slow roll approximation 

First, let us consider the vector sector and eliminate non-dynamical variable /?3 from 
the action for the 2-dimensional vector sector. Now, we define canonically normalized 
variables as 



r = 



b \ k v\ T 
V2k ' 



D = fD. 



(4.16) 



Then, using these canonical variables, we obtain the reduced action for physical variables 



S 



vector 



dr)d 3 k 



V'p i 1 ( m " iA in 2 



2 / 



f 2 v' 2 k 2 ' 



k 2 



\D\ 



+ 1 fv_aky U> D * + f *> D 



y/2 a b k 
where k is time dependent and given by 



(k/b)' 
(k/b) 



(YD* + T*D 



\ 



x + b 2 ( v ) y 



(4.17) 



(4.18) 
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which becomes constant in the isotropic limit a = b. In the isotropic limit a = b, f and 
D represent the cross-mode of gravitational waves and vector waves, respectively. The 
second line in the action ( 14.17P describes how both waves are interacting to each other. 

Next, we use the slow roll approximation to simplify the action. To obtain the 
homogeneous background metric, we integrate the following equations 



H 
IP 



H = 3 l6H 



(4.19) 



by assuming t H /eu <C a /a. The resultant expressions are 

a = (-r])- 1 - 6 ", b = (-rj)- 1 -*"- 1 *" . (4.20) 

In this approximation, the universe shows anisotropic power law inflation. We should 
recall, in the second inflationary phase, the variable I is given by 
, c- 1 



(4.21) 

c 

Note that the range (1, oo) for c corresponds to (0, 1) for I. Using the definition of 1Z, 
we obtain 

f 2 v' 2 



3(-r])- 2 Ie H . 
ur 

From Eq. ( 12. 15ft . the background equation for the vector can be found as 
a 2 J 

From this equation, it is easy to deduce the relation 



[_ 2 _ 3e H +T) H - 2Ie H ] 



where T)h is a slow-roll parameter defined by 

- 2^- {2t H - m ) = 2{2t H - m^-vr 1 ■ 



(4.22) 



(4.23) 



(4.24) 



(4.25) 



Of course, t]h is not related to the conformal time rj. Furthermore, we obtain 

r 
/ 

Substituting these results into the action, we obtain the action in the slow roll 
approximation [72] 



-r])- 2 [2 + 9e H - 3r] H + 61 e H ] . 



(4.26) 



S' 



vector 



dr]d 3 k 



+ 



1 



1. =/ l2 lr 
+ -AD 2 + 



-k 2 + {-riY 2 [2 + 3e H + 3Ie H + 3Ie H sin 2 6} 
-k 2 + {-r,)- 2 {2 + 9e H - 3r, H + 6Ie H sin 2 0) 



iri 



+ ^ m * {-r ] )- 1 sme{T'D* + T*'D) - y^^(- v y 2 sm9(TD* + T*D) 



where we have defined 

• n ^y a 
sin 9 = 

kb 



\D\ 2 
,(4.27) 

(4.28) 
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This 9 represents the direction dependence. In the isotropic limit / = 0, the Lagrangian 
for f becomes the familiar one for gravitational waves in a Friedman-Lemaitre universe. 

In a similar way, we can derive the quadratic action for physical variables in the 
2-dimensional scalar sector. Moreover, it is straightforward to deduce the action in the 
slow roll approximation. The resultant action is given by [72] 



S 



scalar 



dr]d 3 k 



l gg + l jj + l h + l <$>g + L 4>J + l jg 



L 



GG 



J J 



+ 2 

1 

+ 2 



where diagonal parts are given by 

\\o'\ 

- frr 

= -\6(j)'\ 



-k 2 + (-T])- 2 {2 + 3e H + 3Ie H + 3Ie H sin 2 o}] \G\' 
-k 2 + (-T])- 2 {2 + 9e H - 3r] H - QIe H sin 2 6} 



1 

+ 2 



2 + 9e H 



-k 2 + (- V )- 2 
and the interaction parts reads 

L 0G 



3r] H 
1-1 



121 ( 241 



sin 2 fl 



L 4>J 



61 



1-1 



^j(-Vy 2 sin 2 6 {G54>* + G*5 
r])- 1 sin 9 (5(f>*'j + 54>' J* 



(4.29) 

(4.30) 
(4.31) 

(4.32) 
(4.33) 



61 



-T])- 2 sin 9 [5(f)* J + 5(f) J*) 



JG 



(—77) sm9 (G* J + G J 



+ 



H 



r])- 2 sin 9 [G* J + GJ* 



Here, we defined canonical variables 

f\K\ 



G = V2bG , J 



k 



J 



(4.34) 



(4.35) 



(4.36) 



Note that G , J and 5cj) represent the gravitational waves, the vector waves, and the 
scalar perturbations, respectively. The above action shows there exist the interaction 
among these variables. We notice the scalar part ( 14.321) contains / without suppression 
by a slow- roll parameter en- Therefore, to obtain the quasi-scale invariant spectrum of 
curvature perturbation, / itself has to be small. 

From the actions (I4.27P and (I4.29p . we see there are two sources of statistical 
anisotropy of fluctuations. First, the statistical anisotropy of fluctuations comes from 
the anisotropic expansion itself. Intuitively, this can be understood from the anisotropic 
effective Hawking temperature H e g/2n, where H eS denotes the effective expansion rate. 
Indeed, the expansion rate in the direction of the background vector is relatively small, 
hence the effective Hawking temperature is low. Then, this direction has less fluctuation 
power compared to the other directions. Thus, the effective Hawking temperature 
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induces the anisotropy in the power spectrum of fluctuations. This effect is encoded in 
(I4.30p . (I4.3ip . and (I4.32p . The other source of the statistical anisotropy of fluctuations 
comes from the couplings (14. 33p . (14.341) and (I4.35f) due to the background vector field. 
The essential structure of couplings can be understood without complicated calculations. 
Take a look at the following term 

V^g^g^fi&F^ . (4.37) 

Here, we should recall the order of magnitude of background quantities 

/V 2 u ^ i 

For example, to obtain the J — G coupling, one of F^ u have to be replaced by the 
background quantity v'. Hence, the coefficients in the J — G coupling should be 
proportional to fv' which is of the order of y/Ien- This explains the strength of the 
coupling in (I4.35p . Similarly, J — 5(f) coupling should be proportional to fyv' because 
we have to take the variation with respect to <fi. Hence, we can estimate its magnitude 
to be \/J. This explains the interaction term (I4.34[) . Finally, the coupling G — 5(f) has 
a magnitude of the order of f^v' 2 which is proportional to I-JTh. This shows a good 
agreement with the coupling (I4.33[) . Thus, we can understand why there is a hierarchy 
among the couplings of the gravitational waves, the vector waves and the scalar field. 



4-3. Statistical Anisotropy 

In this subsection, we will calculate corrections to power spectrum of various variables 
due to the anisotropy. To set the initial conditions, we need to quantize this system 
by promoting canonical variables to operators which satisfy the canonical commutation 
relations. The point is that, with a given wavenumber, the actions (I4.27P and (I4.29p 
reduce to those of independent harmonic oscillators in the subhorizon limit — kr] ^> 1. 
We choose the Bunch-Davis vacuum state |0) by imposing the conditions a a ,k|0) = at 
an initial time rji. Here, a aik is an annihilation operator whose commutation relations 
are given by 

= 5 ab 5 {3 \k - k'), [a a , k , a 6 , k ;] = . (4.38) 

We are interested in the power spectrum of the scalar perturbations 

(01^(77)^(^)10) = P 50 (k)5(k + p) , (4.39) 

and the power spectrum of the cross and plus mode of gravitational waves 

(0|r k (77)r p (77)|0) ee P r (k)5(k + p) , (4.40) 
(0\G k (r))G p (ri)\0) = P G (k)6(k + p) . (4.41) 

We can also calculate the cross correlation between the plus mode of gravitational waves 
and the scalar perturbations 

(0|^ k (77)G p (77)|0) = P 50G (k)5(k + p) . (4.42) 
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We treat the anisotropy perturbatively and estimate its magnitude by using 
perturbation in the interaction picture. In the interaction picture, the expectation 
value for a physical quantity 0{rj) is given by 



{in \ 0{rj) \ in) 







T exp h / Hi{n)dn 
\ Jm 



0( V ) 



T exp 



H I {n)dn' 



0) ,(4.43) 



where \in) is an in vacuum in the interaction picture, T and T denote a time-ordered 
and an anti-time-ordered product and Hi denotes the interaction part of Hamiltonian 
in this picture. This is equivalent to the following 



{in \ 0{rj) \ in) 



N=0 



dn 



N 



V, 



>IN 



dn 



N-l 



Vi 



'12 



x (0 [Hj{ m ), ■ ■ ■ [Hj{n N ), 0{rj)} ■ • -]]| 0) . (4.44) 

In our analysis, we assume the noninteracting part of Hamiltonian to be that of free 
fields in deSitter spacetime 



2("»?r 2 ) \Qn\ 



and the operators in the interaction picture are given by 
Qn,v{v) = u{r])a n;k + u{n)*a ] n . 



-k- 



u(n) 



-ikr) 



1 kr l/ 



(4.45) 

(4.46) 
(4.47) 



where Q n represent the physical variables D, T, G, J, 5(f). And the rest of the Lagrangian 
f l4.27p -f l4~35]) is regarded as the interaction part Lj = L^> — L . To see the leading effect 
on the anisotropy in the scalar perturbation, which is of the order of J, we evaluate the 
correction due to the interaction given by 



Hf = 



d 3 k 
d 3 k 



Jd 3 k -J-^j{-n)- 1 sm6( y 6(f) r J + 



6/ 



ny 2 sin e (S4) 1 j + 54>j ] 



(4.48) 



Note that in the analogy with the slow-roll parameter in the ordinary slow-roll inflation, 
the term proportional to I sin 2 Odcfidft in ( I4.32p can be expected to give the anisotropy 
8 {in 8(f>k84> p in) / {0 £0k<50 p 0) ~ sin 2 6IN{k) where N{k) is the e-folding number from 
the horizon exit. Thus, the leading correction comes from the interaction through the 
term Hf J . The leading correction is given by 

8 {in 84>k(n)84> p (n) in) 



d V2 I d Vl (0 Hf J ( Vl ), Hf J (ri 2 ),6Mv)8<f> P (v) 



<t>J( 







(4.49) 
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Using Eqs. (14.46)) and commutation relations for the creation and annihilation operators, 
we obtain the anisotropy expressed as follows 

(v) 



5 {in 




in) 


(0 




0) 



24/ 
1-1 
x Im 



sin 9 / dr]2 / dn 



\u(rj)\ 



-Im 



rj 2 ) u (rj 2 )u*(n) + {-n 2 ) u(r) 2 )u*(r)) 



u{v x )u*{r] 2 ) {-(-^i) 1 u(r, 1 )u*\n) + (-r ll ) 2 u(n 1 )u*(n)}] , (4.50) 

where Im denotes the imaginary part. Substituting the function form of u ( I4.47[) and 
introducing time variables x = krj, xi = kn\ an d X2 = ^72, we have 



<5(m 




in) 


(0 




0) 



(x) 



6/ 



sin 2 9 / c?X2 
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XX1X2 

The contribution to the integral from the subhorizon — Xi 3> 1 is negligible. In the limit 
of superhorizon — Xi C 1, we also have — X2 <C 1> — X ^ 1- Hence, the integrand in 
Eq.f l4.5ip approximately becomes 8/xiX2- Thus, the anisotropy can be evaluated as 
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(4.52) 



where N(k) = — ln(—kn) is the e-folding number from the horizon exit. 

For the anisotropy in both polarizations of gravitational waves, the similar 
calculations give 
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(4.53) 



where we used the interaction term in the action (I4.27P for V and that in (14.35)) for G. 
It is interesting to calculate the cross correlation. The leading contribution comes from 
Hj G and Hf . The result is as follows [72] : 

(m|£0 k G p |m) nAT / €h 



(O|5<M0 P |O> 



-24/, 



1 



:N 2 (k) . 



(4.54) 



As we will soon see, this might give a detectable number. 
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5. How To Test Anisotrpic Inflation 

Now, we are in a position to discuss cosmological implication of an anisotropic 
inflationary scenario. As we have listed up in the introduction, there are many 
interesting phenomenology in anisotropic inflation. Here, we recapitulate the results. 
Remember that the anisotropy in the power spectrum is parameterized by 

P(k) = P(k) [l + #*sin 2 fl] . (5.1) 

Then, we can predict the following: 

• There exists statistical anisotropy in curvature perturbations of the order of 

g s = 2AIN 2 (k) . (5.2) 

• There exists statistical anisotropy in gravitational waves of the order of 

g t = QIe H N 2 (k) . (5.3) 

• These exists the cross correlation between scalar perturbations and gravitational 
waves of the order of — 24/^/e/7iV 2 (A;). Using the definition of curvature 
perturbations ( = 5(p/ -*/2e#, one can translate the cross correlation (14.541) between 
the scalar perturbations and gravitational waves to that between the curvature 
perturbations and gravitational waves normalized by the power spectrum of 
curvature perturbations: 

'•-*^-^™-<*»- 

Here, we should note that I <C 1. Due to the interaction on superhorizon scales, there is 
an enhancement factor N 2 (k) in the above quantities. Because of this enhancement, even 
when the anisotropy of the spacetime is quite small, say H/H ~ 10~ 7 in our example, 
the statistical anisotropy imprinted in primordial fluctuations can not be negligible in 
precision cosmology. 

It is useful to notice that there exist consistency relations between observables 

4#i = t H g S) r c = -\/2e H g s . (5.5) 

The consistency relations allows us to test anisotropic inflation in a model independent 
way. Let us explain how to use it. It is known that the current observational limit 
of the statistical anisotropy for the curvature perturbations is given by g s < 0.3 [35J. 
Now, suppose that we detected g s = 0.3. Then, the consistency relations would give 
us predictions. Namely, anisotropic inflation implies the anisotropy in the gravitational 
waves 

9t - 1(T 3 (5.6) 

and the cross correlation 

r c = -V2g s e H ~-4x 10" 3 , (5.7) 

where we used g s ~ 0.3 and e# ~ 10~ 2 . If these predictions are confirmed by the CMB 
observations, that must be a strong evidence of anisotropic inflation. 
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Figure 4. The TT spectra CjJ +2 induced by anisotropy in scalar perturbations, 
that in tensor perturbations and cross correlation. The parameters are chosen as 
a* = 0.3, r = 0.3. 



Apart from the above uses, we can use the observational upper bound g s < 0.3 
to give a cosmological constraint on the gauge kinetic function. Indeed, we have a 
constraint 24IN 2 (k) < 0.3. Since I is derived from the gauge kinetic function and 
the e-folding number N(k) can be determined once reheating process is clarified, the 
constraint on g s implies the constraint on the gauge kinetic function. 

In [71], it is pointed out that the sign of g s predicted by our models is different 
from the observed one. However, it is possible to modify the model so that the sign 
of g s is flipped. For example, we can consider two vector fields. Remarkably, the 
dynamics of vector fields tends to minimize the anisotropy in the expansion of the 
universe and leading to the orthogonal dyad [60] . Then, the orthogonal direction to the 
plane determined by two vectors becomes a preferred direction. In this case, we can 
expect the sign of g s becomes opposite. We can also utilize anti-symmetric tensor fields 
to achieve the same aim. 

We showed in the above how the consistency relations are used to predict the 
observables. So, the next question is how can we see these features in the CMB. 
The answer is that the anisotropy related to tensor perturbations induces off-diagonal 
TB,EB spectra Cg^+i as well as on- and off-diagonal TT, EE, BB,TE spectra C^, 
Ce/ + 2- Here, we have defined the angular power spectrum Ce/> =< a>i, m a t' ,m' > with 
coefficients ae m of the spherical harmonic expansion. 
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Figure 5. The TB and EB spectra Cjf +l , Cfjr +1 induced by the cross correlation. 
As a reference, the conventional diagonal BB spectrum induced by isotropic part of 
the tensor perturbations is plotted with a dotted line. The parameters are taken as 
g s = 0.3, r = 0.3. 



First, we compare the amplitudes of signals induced by the three kind of 
anisotropy [73] . In Fig 0], we have depicted off-diagonal TT correlations CjJ +2 induced 
by the anisotropy. As for the parameter of anisotropy in scalar perturbations, we 
adopted the value g s = 0.3 as a reference, which is just of the order of a systematic 
error in WMAP data. Note that, according to [35J, a signal as small as 2% can be 
detected with the PLANCK. We also assumed the tensor-to-scalar ratio to be r = 0.3. 
Then the other quantities can be determined by the consistency relations in our model: 
r = 16e#, Ag t = e H g s , r c = —\f2e H g s . We see that the contributions of the anisotropy 
in tensor perturbations and the cross correlation are suppressed in comparison to that of 
the anisotropy in scalar perturbations. And, the cross correlation has the contribution 
next to that of the anisotropy in scalar perturbations. This reflects the hierarchy among 
g s , rg t = 0(g s e%), \frr c = 0(g s en)- It is also true for EE and TE spectra. The ratio 
between these effects are given by the slow-roll parameter e# (or the tensor-to-scalar 
ratio r) and does not depend on the value of g s . 

Next we consider peculiar signals of anisotropic inflation [73]. In Fig El we have 
depicted examples of TB and EB correlations Cj^ +l1 Cff +V The parameters are 
again r = 0.3 and g s = 0.3. As a reference, the conventional BB diagonal spectrum 
induced by the isotropic part of tensor perturbations is also plotted with a dotted 
line. Note that unlike parity violating cases for which odd parity correlations Cj B , C§ B 
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exist [HI ESI [75J EZ]) our model predicts even parity correlations such as Cf^ as the 
result of parity symmetry of the system. The ratio of TB correlation induced by cross 
correlation to the isotropic BB correlation is not dependent on tu (or r) for a fixed 
value of g s in our anisotropic inflation model. For the optimistic value of g s ~ 0.3, 
both amplitudes become comparable. This simple order estimation implies that the 
TB signal could be comparable to that of B mode correlation induced by primordial 
gravitational waves. Hence, anisotropic inflation can be a potential source of off-diagonal 
TB correlation, in addition to other effects such as gravitational lensing and our peculiar 
velocity. Since the current constraints on the TB/TE ratio is of the order of 10 -2 [I], 
we need to improve the accuracy by one more order, which might be achieved by the 
PLANCK. 

6. Conclusion 

In this review article, we tried to explain how anisotropic inflation naturally appears in 
supergravity models and how to test anisotropic inflationary models by observations. 
There are several predictions specific to anisotropic inflation. There exists the statistical 
anisotropy in primordial curvature and tensor perturbations. Furthermore, there exists a 
cross correlation between the curvature and tensor perturbations which can be regarded 
as a kind of the statistical anisotropy. Most importantly, there are consistency relations 
between observables in anisotropic inflationary models. This finding gives rise to a 
model independent test of anisotropic inflation. On the other hand, each observable 
gives the information of the specific model. Actually, we have already given the first 
cosmological constraint on the the gauge kinetic function. 

The existence of anisotropic inflation can be regarded as a counter example to 
the cosmic no-hair conjecture. Indeed, in the presence of the gauge kinetic function, 
there could be vector-hair which leads to the anisotropy in the cosmic expansion during 
inflation. Hence, the cosmic no-hair conjecture should be modified appropriately. 
Recently, we have examined inflation with multi-vector fields and found that inflation 
tends to minimize their hair [60] . 

There are several directions to be explored. It is possible to extend the Bianchi 
type I model to other Bianchi type models. Theoretically, it is also interesting to embed 
anisotropic inflation into string theory. Observationally, we need to check the data 
in the CMB more seriously. It is worth investigating how to quantify the statistical 
anisotropy [THJ, [79J [801 EH [82], ]. In fact, our main message is that the statistical 
anisotropy is logically allowed, natural theoretical models exist, and hence observational 
check of the statistical anisotropy needs to be performed. Indeed, there is no apparent 
reason to respect the non-Gaussianity than the statistical anisotropy. 
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